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Using a simple model for nonlinear Landau-Zener tunneling between two energy bands of a Bose- 
Einstein condensate in a periodic potential, we find that the tunneling rates for the two directions 
of tunneling are not the same. Tunneling from the ground state to the excited state is enhanced by 
the nonlinearity, whereas in the opposite direction it is suppressed. These findings are confirmed by 
numerical simulations of the condensate dynamics. Measuring the tunneling rates for a condensate 
of rubidium atoms in an optical lattice, we have found experimental evidence for this asymmetry. 



PACS numbers: PACS number(s): 03.65.Xp, 03.75.Lm 

The phenomenon of Landau-Zener (LZ) tunneling Q 
is a basic quantum mechanical process. It is based on 
the solution of the Schrodinger equation for a two-level 
dynamics when a parameter of the Hamiltonian system 
is time dependent. If at time t — — oo the system is 
prepared in one adiabatic state of the Hamiltonian, the 
time dependence of the Hamiltonian implies that at time 
t = +oo there is a finite probability that the system will 
occupy the other adiabatic state. As far as this tun- 
neling behaviour is concerned, complete symmetry exists 
between the adiabatic states. Variations of the LZ model 
have been studied 0, 0, 3, and an observation of LZ 
dynamics in classical optical systems has been reported 
[5(. More recently, LZ tunneling within a periodic po- 
tential was studied for a nonlinear two-level system in 
which the level energies depend on the occupation of the 
levels 0, • It was discovered that a nonlinearity with a 
positive sign enhances the tunneling probability between 
the ground band and the first excited band. Moreover, 
Niu and coworkers discovered a nonzero LZ tunneling 
probability even in the fully adiabatic limit when the 
nonlinearity was larger than a critical value pj. Critical 
values for deformations of the energy level structures were 
obtained in refs. In a Bose-Einstcin condensate 

inside a periodic potential such as an optical lattice, the 
mean-field interaction between the atoms can be compa- 
rable to other energy scales of the system and hence the 
level-dependent energy shift can lead to an observable 
modification of the tunneling behaviour. 



In the present work we explore, theoretically and ex- 
perimentally, the Landau-Zener tunneling between Bloch 
bands of a Bose-Einstein condensate in an accelerated op- 
tical lattice. The optical lattice depth controls the tun- 
neling barrier, while the optical lattice acceleration con- 
trols the time dependence of the Hamiltonian. We show 
that the mean-field nonlinearity produces an asymmetry 
for the tunneling probability between the adiabatic states 
of the Hamiltonian. More precisely the tunneling proba- 
bility from the lower energy adiabatic state to the upper 
one is enhanced, while the tunneling probability from the 
higher energy level to a lower level is suppressed. Numer- 
ical integration of the one-dimensional Gross-Pitaevskii 
equation and a simple two-state model demonstrate this 
asymmetry. Moreover, our experimental data obtained 



with a rubidium Bose-Einstcin condensate confirms this 
prediction. 

The qualitative explanation of this asymmetry in the 
tunneling transition probabilities depends on the choice 
of the quantum representation used for the description 
of the nonlinear system. We have identified two differ- 
ent basic mechanisms, acting simultaneously or indepen- 
dently. In the first mechanism the nonlinear term of 
the Schrodinger equation acts as a perturbation whose 
strength is proportional to the energy level occupation. 
If the initial state of the condensate in the lattice cor- 
responds to a filled lower level of the state model, then 
the lower level is shifted upward in energy while the up- 
per level is left unaffected. This reduces the energy gap 
between the lower and upper level and enhances the tun- 
neling. On the contrary, if all atoms fill the upper level 
then the energy of the upper level is increased while the 
lower level remains unaffected. This enhances the energy 
gap and reduces the tunneling. The overall balance leads 
to an asymmetry between the two tunneling processes. 
In a different representation the nonlinearity term in the 
Schrodinger equation produces an additional term to the 
optical lattice depth. This additional term increases or 
decreases the optical lattice depth depending on whether 
the tunneling process proceeds from the lower level to 
the upper one, or vice versa. Thus the tunneling due 
to a weaker or stronger potential barrier is increased or 
decreased, as the case may be, depending on the initial 
state. 

The motion of a Bose-Einstein condensate in an ac- 
celerated ID optical lattice is described by the Gross- 
Pitaevskii equation 
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where M is the atomic mass, Ul — n/d is the optical 
lattice wavenumber with d the optical lattice step, and 
Vq is the strength of the periodic potential depth. The 
s-wave scattering length a s determines the nonlinearity 
of the system. Equation ^ is written in the comoving 
frame of the lattice, so the inertial force Ma^ appears as 
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a momentum modification. The wavefunction ip is nor- 
malized to the total number of atoms in the condensate 
and we define no as the average uniform atomic density. 
Defining the dimensionless quantities E rec = fi 2 k\/2M , 
x = 2k^x, t — 8E rec t/h, and rewriting ip = tp/^Jno, 
v = V /16E rec , a = Ma L /16E rec k L , C = 7ra s n /fcf,, Eq. 
(JIJ is cast in the following form 



comes from the linear effect. We confirmed the presence 
of a tunneling asymmetry by integrating directly Eq. i[TJl 
(taking into account the full experimental protocol de- 
scribed below), finding qualitative agreement with the 
prediction of the two-state model. For small C values, 
we have fitted the C dependence of the tunneling rate of 
Figffl (a) through the following expression 0, Q : 



.dtp If d 



ip + vcos(x)ip + C\ip\\p (2) 



where we have replaced i with x, etc. In the neigh- 
borhood of the Brillouin zone edge we can approximate 
the wave function by a superposition of two plane waves 
(the two level model of ref. 0), assuming that only 
the ground state and the first excited state are popu- 
lated. We then substitute ip(x, t) = a{t)e lqx +b(t)e 



i(q—l)x 



with \a(t)\ 2 + \b(t)\ 2 = 1 in Eq. 
the coefficients of e lqx and e 1 ' 9 " 1 ^, 
netic terms and dropping the irrelevant constant energy 
1/8 + C[l + (H 2 + |6| 2 )/2], Eq. © assumes the form 
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where cr, i = 1, 2, 3 are the Pauli matrices 10] . The adi- 
abatic energies of Eq. |3]have a butterfly structure at the 
band edge of the Brillouin zone for C > v 0, @, , but 
in the present work we always work in a regime where 
C <C v, hence that structure plays no role. 

In the linear regime (C = 0), evaluating the transition 
probability in the adiabatic approximation, we find the 
linear LZ formula for the tunneling probability r 

r = e 2^ (4) 

expressing the occupation changes in terms of the rate a 
at which the diagonal energies of the linear Hamiltonian 
change their value, and of the off-diagonal interaction en- 
ergy v. In the nonlinear regime, as the nonlinear param- 
eter C grows, the lower to upper tunneling probability 
grows as well until an adiabaticity breakdown occurs at 
C = v 7}. The upper to lower tunneling probability, 
on the other hand, decreases with increasing nonlincar- 
ity |ll| . We derived the tunneling rate from the numeri- 
cal integration of Eq. In Fig. ^a) we plot the lower 
to upper tunneling rate (initial (a, b) =(1, 0)) and the up- 
per to lower tunneling rate (initial (a, 6) =(0, 1)) of the 
Bose-Einstein condensate as a function of the nonlinear 
parameter C for different accelerations of the optical lat- 
tice. We see that for C = the rate is the same for both 
tunneling directions whereas for C ^ the two rates are 
different, and the smaller the acceleration the larger the 
difference. This result is intuitive since for very small 
accelerations the main contribution originates from the 
nonlinear effect (the linear tunneling of Eq. being 
small) while for large accelerations the main contribution 
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with (3 = 0.75, 0.17, 0.14 for the different acceleration val- 
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FIG. 1: LZ tunneling rate r within the two-level model as a 
function of the nonlinear parameter C for different accelera- 
tions. The rates coincide for C — 0, whereas for C / they 
differ significantly. Results are for v = 0.134 corresponding to 
Vo = 2.2E rec , in (a) for a = 0.01, 0.04, 0.07, corresponding 
to aj, = 0.8, 3.2, 5.6 ms -2 from bottom to top and in (b) for 
a — 0.036 corresponding to a_L = 2.9 ms -2 . In (b) experi- 
mental results are denoted by open symbols for the excited 
band to ground state band tunneling and by filled symbols for 
the ground state band to excited band tunneling. The exper- 
imental points at C = 0.025 have been re-scaled from Fig. [2] 
using the Landau-Zener formula to match the acceleration to 
the one used to obtain the other experimental points. The 
shaded area represents the confidence region for the predic- 
tion of the two-level model taking into account the uncertainty 
in our measurement of the lattice depth. 



The nonlinear regime is more easily interpreted by 
writing Eq. (©) as 
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Here the nonlinear diagonal terms represent a shift of the 
energy levels, different for an occupation of the ground 
or upper level, decreasing or increasing the energy gap 
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depending on the initial state. The off-diagonal terms 
of Eq. produce an equivalent contribution to the non- 
linear tunneling process. The off-diagonal terms modify 
the interaction term v equivalent to a Rabi frequency 
in the two-level model. The off-diagonal scalar product 
between the two states a* b, calculated applying the adi- 
abatic approximation technique of |12| . changes sign de- 
pending on the initial state (lower or upper). Thus the 
linear Rabi frequency is modified by the nonlinear off- 
diagonal term, and the increase or decrease depends on 
the initial state. This nonlinear modification of the in- 
teraction energy constitutes an additional contribution 
to the asymmetric tunneling rate. 

Experimentally, we investigated the phenomenon of 
asymmetric tunneling between the energy bands of Bose- 
Einstein condensates in an optical lattice using a setup 
described in detail in 

[HH|. Briefly, we create con- 
densates of N ps 1 4 rubidium atoms in a time-orbiting 
potential (TOP) trap. Once condensation has been 
achieved, the mean trapping frequency v tr ap of the mag- 
netic trap is adiabatically reduced to values between 
15 Hz and 50 Hz. Thereafter, two laser beams with waists 
of 1.8 mm and intersecting at an angle 9 = 38deg at the 
position of the condensate are switched on with a linear 
ramp of duration T ramp = 10 ms, thus ensuring adiabatic- 
ity of the loading process. The beams are detuned to the 
red side of the rubidium atomic resonance by ps 30 GHz 
and have a variable frequency difference Av between 
them, controllable through two acousto-optic modulators 
which are also used to vary the intensity of the beams. 
In this way, a periodic potential with lattice constant 
d = 1.18 /im and lattice recoil energy E rec /h = 455 Hz is 
created, which through the frequency difference Ais can 
be made to move at a constant velocity v — dAv or ac- 
celerated with cll = d^f-- For the current experiment, 
lattice depths between 0.25 E rec and 2.5 E rec were used. 

Landau-Zener tunneling between the two lowest energy 
bands of a condensate inside an optical lattice is inves- 
tigated in the following way. Initially, the condensate is 
loaded adiabatically into one of the two bands. Subse- 
quently, the lattice is accelerated in such a way that the 
condensate crosses the edge of the Brillouin zone once, 
resulting in a finite probability for tunneling into the 
other band (higher-lying bands can be safely neglected as 
their energy separation at the edge of the Brillouin zone 
is much larger than the band gap). After the tunnel- 
ing event, the two bands have populations reflecting the 
Landau-Zener tunneling rate (assuming that, initially, 
the condensate populated one band exclusively). In or- 
der to experimentally determine the number of atoms 
in the two bands, we then increase the lattice depth 
(from ps 2 E rec to ps 4 E rec ) and decrease the accelera- 
tion (from ?s 3ms" 2 to ps 2ms -2 ). In this way, suc- 
cessive crossings of the band edge will result in a much 
reduced Landau-Zener tunneling probability between the 
ground state band and the first excited band (of order a 
few percent). The fraction of the condensate that af- 
ter the first tunneling event populated the ground state 



band will, therefore, remain in that band, whereas the 
population of the first excited band will undergo tunnel- 
ing to the second excited band with a large probability 
(around 90 percent) as the gap between these two bands 
is smaller than the gap between the two lowest bands by 
a factor ps 5 for our parameters. Once the atoms have 
tunneled into the second excited band, they essentially 
behave as free particles since higher-lying band-gaps are 
smaller still, meaning that the fraction of the conden- 
sate that populated the first excited band after the first 
tunneling event will no longer be accelerated by the lat- 
tice. In summary, using this experimental sequence we 
selectively accelerate that part of the condensate further 
that populates the ground state band. In practice, in or- 
der to get a good separation between the two condensate 
parts after a time-of-flight, we accelerate the lattice to a 
final velocity of 4 — 6 v rec [15j and absorptively image the 
condensate after 22 ms (see Fig- EJ • 
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FIG. 2: Profiles of absorption images taken after 22 ms of 
time-of-flight of condensates released after the acceleration 
procedure described in the text. The condensates were pre- 
pared in the ground state band in (a) and in the first excited 
band in (b) within an optical lattice with depth Vb = 2.6 E rec 
and ol = 2.9 ms -2 . 

In order to investigate tunneling from the ground state 
band to the first excited band, we adiabatically ramped 
up the lattice depth with the lattice at rest and then 
started the acceleration sequence. The tunneling from 
the first excited to the ground-state band is investigated 
in a similar way, except that in this case we initially pre- 
pare the condensate in the first excited band by moving 
the lattice with a velocity of 1.5v re c (through the fre- 
quency difference Av between the acousto-optic modula- 
tors) when switching it on. In this way, in order to con- 
serve energy and momentum the condensate must popu- 
late the first excited band at a quasi-momentum half-way 
between zero and the edge of the first Brillouin zone [lj| • 
Thereafter, the same acceleration sequence as described 
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above is used in order first to induce Landau-Zener tun- 
neling and then to separate the fraction of the conden- 
sate that underwent tunneling to the lowest band from 
the one that remained in the first excited band. For both 
tunneling directions, the tunneling rate is measured as 



N, 



tunnel 

N tnt ' 



(7) 



where N to t is the total number of atoms measured from 
the absorption picture. For the tunneling from the first 
excited band to the ground band, N tunne i is the number 
of atoms accelerated by the lattice, i.e. those detected in 
the final velocity class 4iv e c, whereas for the inverse tun- 
neling direction, Ntunnei is the number of atoms detected 
in the v = velocity class. To illustrate our method 
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FIG. 3: Landau-Zener tunneling between the two lowest en- 
ergy bands of a condensate in an optical lattice as a function 
of the lattice depth (carrying a ±10% systematic error) for 
a = 0.025 corresponding to «i = 2.0 ms -2 . Tunneling rates 
from the ground state band to the first excited band (filled 
symbols) and vice versa (open symbols) are virtually identical 
and agree with the linear prediction (dashed line) as the non- 
linear parameter C ~ 0.025 is small for the trap used in these 
measurements (77 trap = 15 Hz). The deviation from theory of 
the filled symbol for small lattice depths is most likely due 
to a breakdown of the adiabatic approximation of the linear 
Landau-Zener model. 



for measuring the tunneling in both directions, in Fig. [3| 
we report the two tunneling rates as a function of lat- 
tice depth for a condensate in a weak magnetic trap and 
hence a small value of the interaction parameter C ^3 • 



In this case, both tunneling rates are essentially the same 
and agree well with the linear Landau-Zener prediction. 
By contrast, when C is increased, the two tunneling rates 
begin to differ, as can be seen in Fig.^b). Qualitatively 
we find agreement with the theoretical predictions of the 
non-linear Landau-Zener model, whereas quantitatively 
there are significant deviations. We believe these to be 
partly due to experimental imperfections. In particular, 
the sloshing (dipolar oscillations) of the condensate inside 
the magnetic trap can lead to the condensate not being 
prepared purely in one band due to non-adiabatic mixing 
of the bands if the initial quasimomentum is too close to 
a band-gap. Furthermore, a numerical simulation of the 
experiment shows that for large values of C, for which the 
magnetic trap frequency was large, the measured tunnel- 
ing rates are significantly modified by the presence of the 
trap. We have, however, verified that when C in the sim- 
ulation is varied without varying the trap frequency, the 
asymmetric tunneling effect persists. 

In summary, we have numerically simulated Landau- 
Zener tunneling between two energy bands in a periodic 
potential and found that, in the presence of a nonlinear 
interaction term, an asymmetry in the tunneling rates 
arises. Experimentally, we have measured these tunnel- 
ing rates for different values of the interaction parameter 
and found qualitative agreement with the simulations. 
We note here that the phenomenon of asymmetric tun- 
neling should be a rather general feature of quantum sys- 
tems exhibiting a nonlinearity. For instance, calculating 
the energy shift due to a nonlinearity for two adjacent lev- 
els of a harmonic oscillator, one finds that both levels are 
shifted upwards in energy, the shift being proportional to 
the population of the respective level. The energy differ- 
ence between the levels, therefore, decreases if only the 
lower state is populated and increases if all the popula- 
tion is in the upper level. Finally, we note that state- 
dependent mean-field shifts have also been observed in 
measurements of the clock-shift in ultra-cold and Bose- 
condensed atomic samples fl9| . 
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